Spectral analysis of boundary layers in Rayleigh-Benard convection 
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A combined experimental and numerical study of the boundary layer in a 4:1 aspect-ratio 
Rayleigh-Benard cell over a four-decade range of Rayleigh numbers has been undertaken aimed 
at gaining a better insight into the character of the boundary layers. The experiments involved the 
simultaneous Laser Doppler Anemometry (LDA) measurements of fluid velocity at two locations, 
i.e. in the boundary layer and far away from it in the bulk, for Rayleigh numbers varying between 
1.6 X 10^ and 2.4 x 10^. In parallel, direct numerical simulations (DNS) have been performed for 
the same configuration for Rayleigh numbers between 7.0 x 10** and 7.7 x 10^. The temperature and 
velocity probability density functions and the power spectra of the horizontal velocity fluctuations 
measured in the boundary layer and in the bulk flow are found to be practically identical. Except 
for the smallest Rayleigh numbers, the spectra in the boundary layer and in the bulk central region 
are continuous and have a wide range of active scales. This indicates that both the bulk and the 
boundary layers are turbulent in the Ra number range considered. However, molecular effects can 
still be observed and the boundary layer does not behave like a classical shear-driven turbulent 
boundary layer. 

PACS numbers: 47.20.Bp, 47.27.nb, 44.25. +f 

Keywords: Rayleigh-Benard convection, turbulence, boundary layer, DNS, LDA, power spectrum 



I. INTRODUCTION 



Rayleigh-Benard (KB) convection in a fluid trapped 
between two horizontal plates of unequal temperature 
with the bottom wall being warmer than the top wall, 
has long served as a paradigm of thermal convection. 
Despite numerous studies (see e.g. overviews by [H, 01), 
it continues to pose challenges because of a number of 
still unresolved and controversial issues. We recall that 
RB convection is characterised by the Rayleigh number 
Ra = agAQH^ / (i'k), and the Prandtl number Pr = i^/k, 
where a is the thermal expansion coefRcient, AQ is the 
temperature difference between the cold and the hot 
plate, H is the vertical distance between the plates, v 
is the kinematic viscosity, and n is the thermal diffusiv- 
ity. In practical RB set-ups the lateral size, W, of the 
plates is finite. This is taken into account by a third di- 
mensionless number known as the aspect ratio F — W/H. 

A topic of intense debate is the scaling law for the di- 
mensionless heat transfer, Nu ~ Ra"*^, especially for the 
range Ra > 10^ Here Nu is the Nusselt number 
Nu = hH/X = -iJ/(Ae)ae/52;|waii, h stands for the 
convective heat transfer coefficient, A is the thermal con- 
ductivity, and 9 is the averagedtemperature. 0, H, 0| 
proposed scaling laws for the heat transfer and tempera- 
ture statistics in RB convection, which differed from the 
classical theories, such as that of p]. In the search for 
general scaling behaviour, Q presented a unifying the- 
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ory dividing the Ra, Pr plane into several regions, each 
with different scaling. Crucial assumptions in this theory 
are the existence of a large-scale circulation, also known 
as the wind, and a Blasius type laminar boundary layer. 
The assumption that the boundary layers of Rayleigh- 
Benard convection behave as a laminar Blasius boundary 
layer is supported by several observations: (i) the scal- 
ing of the friction factor C/ is compatible with that of 
a laminar boundary layer past a flat plate at moderate 
Ra d, [13, [HI; (ii) the typical Reynolds number Re of 
the flow is low at moderate Ra: for water Re ~ 200 at 
Ra 10^ and F = 4, which is generally accepted to be 
too low to sustain turbulence. p\ argued that at large 
enough Ra this laminar boundary layer breaks down and 
becomes turbulent, presuming that the transition to a 
turbulent boundary layer is shear triggered. The critical 
Ra number for this transition was estimated to be in the 
order of 10^** for Pr = 1. 

Despite this apparently convincing evidence, it is not 
clear how a boundary layer which is highly unsteady 
due to continuous plume impingement and detachment 
can behave quantitatively as a laminar boundary layer. 
This would indicate that the plumes only cause unsteady 
disturbances and a consequent time-dependence in the 
boundary layer but are passive otherwise. 

The aim of this paper is to enhance the understand- 
ing of the dynamics of the boundary layers by comparing 
velocity statistics in the boundary layer with those de- 
termined in the bulk. It has widely been regarded that 
the main prerequisite for unsteady velocity fluctuations 
to be qualified as turbulence, in contrast to other un- 
correlated, chaotic and random, fluctuations, is a power 
spectrum that is continuous over a wide range of scales. 
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For this reason we focused on the analysis of spectral den- 
sity functions (sdf), complemented with the probability 
density functions (pdf), aimed at clarifying the nature of 
the fluctuations in the boundary layer. For this purpose 
we use in parallel complementary numerical and exper- 
imental techniques: direct numerical simulations (DNS) 
are performed for the lower Ra values, Laser Doppler 
Anemometry (LDA) measurements are conducted for the 
medium and higher Ra values. All experiments and sim- 
ulations combined cover a Ra range between 7.0 x 10''' 
and 2.4 x 10^. In the midrange both techniques overlap 
over about one decade. 



II. METHODOLOGY 
A. Experiments 

The experiments were conducted in a 600 x 600 x 
155 mm'^ {W xWxH) cell filled with water (described in 
more detail in [11]). In brief, two copper plates at the top 
and bottom are kept at a constant temperature by pass- 
ing water through the plates' internal channels, drawn 
from two basins containing constant temperature water 
(inaccuracy less than 0.03 K). The plates thus impose 
a controlled temperature difference AQ on the working 
fluid. Throughout this paper we will use a Cartesian 
coordinate system with its origin at the centre of the 
bottom wall, Fig. [TJ The x- and y- coordinates are mea- 
sured along the bottom wall, parallel to the side walls of 
the RB cell. The z-coordinate measures the wall normal 
distance. 

Velocities were measured by using two one-component 
laser Doppler anemometers manufactured by Dantec. 
The green (514.5 nm) and blue (488.0 nm) colours of 
a 5 W Argon-ion laser are used to simultaneously mea- 
sure the horizontal component in the x-direction in the 
bulk and in the boundary layer. The measured velocity 
component is in the direction of the large-scale circu- 
lation. To enable the detection of instantaneous flow 
reversals one beam of each colour was frequency pre- 
shifted by a Bragg cell. The light scattered by 5 ^m 
polyamid seed particles was collected in on-axis forward 
direction, which, in comparison to back-scatter, results 
in a superior signal-to-noise ratio of the Doppler signals. 
The photomultiplier output signals were first downmixed 
electronically to an effective pre-shift of 30 kHz and then 
fed to two Dantec BSA signals processors to determine 
the instantaneous velocities. The LDA system measured 
the instantaneous velocities of the seed particles with an 
inaccuracy of 0.5%. The length and diameter of both 
measurement volumes are 0.8 mm and 0.1 mm, respec- 
tively, which is sufficiently small to resolve the smallest 
length scales of the flow. For each Ra number, the mea- 
surement time T was at least 20 hours to have sufficient 
data for an accurate statistical description of the flow. 

The mean data rate varied between 10 Hz and 30 Hz, 
depending on Ra and distance from the wall. For mod- 
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FIG. 1: Schematic of the measurement points in the RB cell 
for both DNS and LDA measurements 
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FIG. 2: Kinematic boundary-layer height A,,, based on Equa- 
tion [5] as well as from DNS, and the measurement locations, 
Zhi, as function of Ra. 



erate Ra numbers, the relatively low characteristic veloc- 
ities result in a relatively low data rate. Measurement 
locations close to the wall have an inherently low data 
rate due to the low fluid volume flux. 

Table U shows the experimental parameters as well as 
the relevant non-dimensional numbers. The Rayleigh 
number is varied by changing the temperature difference 
A6 between the bottom and top plates from nearly 1 K 
at Ra = 1.6 x 10^ to 25 K at Ra = 2.4 x 10^. For 
the highest Ra numbers, the relatively large tempera- 
ture differences cause larger density differences, but these 
remain small enough to neglect non-Boussinesq effects. 
The Prandtl number will change as well, but the ratio 
between the smallest and largest Pr number is only 1.7. 
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Ra 


Pr 


Ae [K] 


e [°c] 


a [K-i] 




V [m^ s 




T [s] 


1.6 X 10^ 


9.4 


0.92 


10.2 


8.8 X 10" 


-5 


1.3 X 10^ 


-6 


1.1 X 10^ 


8.2 X 10^ 


7.0 


1.57 


19.9 


2.1 X 10" 


-4 


1.0 X 10 


-6 


0.9 X 10^ 


4.2 X 10** 


5.5 


6.52 


30.0 


3.0 X 10" 


-4 


8.0 X 10^ 


-7 


1.2 X 10^ 


2.4 X 10^ 


5.5 


25.1 


30.3 


3.0 X 10" 


-4 


8.0 X 10 


-7 


1.3 X 10^ 



TABLE I: Experimental details for the LDA measurements. All experiments have been performed in a F = 4 aspect ratio 
domain. T denotes the measurement time. 



B. Simulations 

The direct numerical simulations cover the lower to 
moderate Ra number range between Ra = 7.0 x 10'* and 
Ra = 7.7 X 10^. The numerics are based on staggered 
finite differences, with a central second order scheme for 
advection, and a second order Adams-Bashforth time in- 
tegration scheme. To effectively resolve the boundary 
layers, grid stretching in the wall-normal direction has 
been implemented, and special care has been taken to 
preserve the skew-symmetry of the advective operator 
jist - The code, as reported in [l4||, has been extended to 
support sidewalls. In the situation without sidewalls, the 
Poisson equation is Fourier-transformed in the two ho- 
mogeneous directions. As a result, the Poisson equation 
decouples for these directions, and the pressure ampli- 
tudes can be solved per component by solving a tridiag- 
onal system in the wall-normal direction. To incorporate 
sidewalls, the Fourier-transform was replaced by a co- 
sine transform, which automatically satisfies a Neumann 
boundary condition for pressure. The sidewall bound- 
ary conditions for the velocity and temperature could be 
enforced directly due to the explicit time-integration. 

Dirichlet boundary conditions are applied at the top 
and bottom plate, i.e. no-slip and isothermal boundary 
conditions for velocity and temperature, respectively. On 
the sidewalls, Neumann boundary conditions are applied 
for both temperature and velocity, i.e. free-slip and insu- 
lating walls for velocity and temperature. The free-slip 
(stress free) velocity boundary conditions were used to 
avoid having to resolve the boundary layer at the side- 
walls. Although this changes the fiow physics near the 
sidewalls, its global effect is limited for large aspect ra- 
tio domains, as preliminary simulations showed. Neu- 
mann boundary conditions were applied for pressure at 
all boundaries. 

In the simulations, the properties of water are used 
except that Pr = 7 is kept fixed for all Ra. The grid res- 
olution (Table ini is sufficient to resolve the smallest tur- 
bulent scales, i.e. the Kolmogorov scale ?7k = {v'^ /sY^^ 
and the Corrsin scale rjc = Pr~*/^ ?7k- At Pr = 7, the 
Corrsin scale dominates the grid requirements, and the 
resolution roughly corresponds to Aa; ~ Ay « Az « 477c- 
To resolve the boundary layers at the top and bottom 
wall properly, grid stretching has been applied in the 
wall-normal direction, such that at least eight points 



Ra 


grid 


M/t* 




T/t* 


Nu 


7.0 X 10'' 


64^ X 48 


4.0 X 10" 


-4 


240 


3.8 


5.6 X 10^ 


96^ X 64 


2.3 X 10" 


-4 


140 


7.2 


2.6 X 10^ 


128^ X 92 


1.7 X 10" 


-4 


165 


11.3 


7.2 X 10^ 


192^ X 128 


1.3 X 10" 


-4 


155 


15.2 


2.1 X 10^ 


384^ X 160 


1.0 X 10" 


-4 


40 


20.0 


7.7 X 10^ 


384^ X 160 


1.0 X 10" 


-4 


30 


30.0 



TABLE II: Simulation details for DNS. 



are within the thermal boundary layer. It should be 
noted that the simulation at the highest Ra number 
{— 7.7 X 10^) may be slightly underresolved in the hor- 
izontal directions. High quality simulation as reported 
in [l^ at Ra — 10* and Pr = 7 required 768 points 
per horizontal direction, but reduction to 384 points was 
the only way to simulate over a sufficiently long time 
interval to obtain reliable one-point statistics required 
for the comparison with the experiments. The addi- 
tional information shown in Table Hi] are the dimension- 
less timestep t/t*, the number of simulated turnovers 
T/t* and the average Nusselt number Nu. The typical 
convective turnover time t* = H/U with U the typical 
velocity based on the free fall scaling. The average Nus- 
selt number is calculated as Nu = —H/{A<d) dQ/dz\^^^^, 

with Q averaged over the entire bottom surface and over 
many turnovers. 



C. Measurement locations 

Simultaneous measurements have been performed in 
the boundary layer and in the bulk of the cell for various 
Ra. The measurement location in the bulk is at a;buik = 
W/3, j/buik = W/2 with a fixed wall distance Zbuik = H/2, 
see Fig. [1] The measurement location in the boundary 
layer has the same x- and y-coordinates, but the wall 
distance is different, i.e., Zbi = 0.4A„, where A„ is the 
kinematic boundary layer thickness determined as in 0] : 

A„ = H/VR^. (1) 
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In this expression for A^, 
taken from 



the Reynolds number, Re, was 



Re = 0.073Rai/2pj.-5/6^ 



(2) 



Table IIIII lists the physical measurement locations, 
where the DNS values have been rescaled to H = 
155 mm. Figure [5] shows that the wall distance for the 
boundary layer measurements slightly deviates from lin- 
ear behaviour (on a double-log plot) due to minor changes 
in Pr for the experimental points. In separate experi- 
ments, LDA was used to determine mean-velocity and 
rms- velocity profiles above the wall (at Ra « 8.0 x 10*). 
The results of these experiments confirmed that the mea- 
surement locations at z = Zbi were all located in the 
boundary layer for the full Ra number range considered 
in this study. The DNS values for A„, denoted by solid 
circles in Fig. [21 support this observation. Clearly, ([T|) 
overestimates the scaling exponent of A„, so that the 
prefactor a in the BL measurement points z^jj — aA„ is 
a decreasing function of Ra. Hence, as Ra increases, the 
measurements are performed deeper into the kinematic 
boundary layer. The disparity between the theoretical 
prediction ([T]) and the measurements for A„ will be dis- 
cussed in more detail at the end of section IIII Bl 

The DNS time series are short in respect to the LDA 
series. To enable a comparison between LDA and DNS 
results, the DNS data are extracted from four points 
in the domain. These points are symmetrically placed 
around the origin with horizontal coordinates {ztW/6, 0), 
(0, ±.W/Q), as shown in Fig.[T] The four obtained time se- 
ries are regarded as nearly statistically independent, and 
are then averaged to effectively improve the statistics of 
the DNS time series. 



III. RESULTS 
A. Probability density functions 

Figure [3] shows the measured and computed variance 
of the horizontal velocity ct^ in the boundary layer and in 
the bulk. The graph shows that for increasing Ra there 
is a gradual transition of the variances between DNS and 
LDA results giving confidence in the procedure followed. 

Temperature pdfs obtained from DNS for the bound- 
ary layer and the bulk are depicted in Figure |4(b)| and 
show an exponential distribution even for Ra = 5.6 x 10^. 
The straight lines (on logarithmic scale) are generally as- 
sociated with the hard turbulent regime, which sets in at 
Ra = 4 X lO'^ for aspect-ratio- 1 cells (e.g. 0, H, [l7j). 
Recall that the transition between soft and hard tur- 
bulence takes place at relatively low Ra number in the 
aspect-ratio-4 cell, see [iSj. The boundary layer pdf is 
skewed towards positive temperatures. This is generally 
attributed to the interaction of thermal plumes with the 
flow. Plumes leaving the boundary layer result in positive 
temperature excursions. 
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FIG. 3: Variances of the horizontal velocity component, u, in 
boundary layer (a) and in the bulk (b) 



Figure [S] depicts the probability density functions 
(pdfs) of the computed and measured u component. As 
expected, the pdfs of both horizontal velocity compo- 
nents were found to have the same shape in the simu- 
lations. For brevity only the u component is therefore 
depicted in Figure [5] For a better comparison, the pdfs 
are scaled with the rms value (t„. The results show that 
the pdfs in both the boundary layer and in the bulk ex- 
hibit a nearly Gaussian-like shape for the complete Ra 
number range. Only the pdfs at Ra = 7 x 10'* deviate 
largely from this trend, mainly due to the presence of a 
strong oscillatory component in the velocity time series, 
as will be shown later when discussing the power spectra. 

[ig. I20I make a distinction between sub-Gaussian, 
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Ra 


Pr 


Re 


Xu/H 






7.0 X 10-* 


7.0 


3.8 


0.512 


0.204 


DNS 


5.6 X 10^ 


7.0 


11 


0.305 


0.121 


2.6 X 10® 
7.2 X lO** 


7.0 
7.0 


23 
39 


0.207 
0.161 


0.0827 
0.0641 




2.1 X 10^ 


7.0 


66 


0.123 


0.0490 




7.7 X 10'^ 


7.0 


1.3 X 10^ 


0.0884 


0.0354 


DA 


1.6 X 10'^ 


9.4 


48 


0.145 


0.0578 


8.2 X lO'^ 


7.0 


1.2 X 10^ 


0.0897 


0.0357 




4.2 X 10** 


5.5 


4.2 X 10^ 


0.0485 


0.0193 




2.4 X lO'' 


5.5 


1.0 X 10^ 


0.0309 


0.0123 



TABLE III: Measurement heights for the different Ra experiments. The Re number is based on equation O 
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FIG. 4: Probability den sity functions for temperature in 
the boundary layer (a) and in the bulk |(b)| of the RB cell 
obtained from DNS. Histograms are scaled with their rms 
value (Te- 



Gaussian and hyper-Gaussian stages in the development 
of turbulence. They argue that the nearly Gaussian 
shape of the pdfs in the present study indicates that 
turbulence is in a fully developed stage. The pdfs for 
the highest Ra values have non-zero mean values. This 
is caused by the large-scale circulation or "wind" within 
the RB cell. For Ra — 4.2 x 10® the mean value is approx- 



imately 5 mms~^. The wind is rather weak in the sense 
that its mean value is small compared to the magnitude 
of the rms velocity fluctuations^ in accordance with other 
findings reported in [l^, Uli (22|, . It is also observed 
that the wind is not always in the same direction. In 
some experiments (Fi g. [31[h)[) it appears to be clockwise, 
while in others (Figs, [^fij y)) it is counterclockwise. 



B. Spectral density functions 

The obtained velocity time series are investigated fur- 
ther by computing the spectral density function (sdf). 
The sdf is defined as the Fourier transform of the au- 
tocorrelation function p(t) = u' {t)u' {t + t) / u'^ , where 
u'{t) (= u{t) — u) is the fluctuating component of the 
turbulent velocity, t is time, r is the time lag, and the 
overline denotes ensemble averaging. All sdfs are cal- 
culated using the same algorithm, to make a fair com- 
parison for both DNS and LDA data. This algorithm, 
described in [23|, can account for the randomness of the 
sampling times in LDA measurements. 

Figure [S] shows the sdfs of the horizontal component 
in the a;-direction determined in the boundary layer and 
in the bulk for the full range of Ra considered in this 
study. The sdfs determined from DNS data are accurate 
at high frequencies, but less reliable in the low frequency 
range due to the relatively small measurement times. In 
contrast, the low frequency region in the LDA sdfs is 
accurate, while the high frequency range suffers from a 
relatively high statistical scatter caused by the random 
sampling process in LDA. The level of the statistical scat- 
ter (which strongly depends on the mean data rate) sets a 
lower limit to the spectral density that can be accurately 
determined from LDA data. 

When comparing the various sdfs in Figure [Slit is clear 
that the sdf for Ra = 7 x 10"* is different from all others. 
At this low Ra numbers, the flow is not turbulent but in a 
state of spatio-temporal chaos. As mentioned earlier, the 
velocity time series is characterised by strong oscillations. 
As a consequence, the sdf in Fig. 6(a) (and to a lesser 
extent the sdf in Fig. |6(b)[ ) is dominated by a limited 
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FIG. 5: Probability density functions of the u horizontal velocity component in the boundary layer (bl) and in the bulk of the 
RB cell. Histograms are scaled with their rms value a-u- The source of the data, DNS or LDA, is indicated. 



number of discrete peaks and the sdf rolls off rapidly at 
higher frequencies. 

Peaks in the sdf also occur in the very low frequency 

The peak frequencies cjp 



range in Figs. 6(i) and 6(j) 



correspond to periods of 27r/tjp = 283 s and 172 s at 
Ra = 4.2 X 10^ and Ra = 2.4 x 10^, respectively. These 
oscillations are unrelated to the ones observed at low Ra, 
and have their origin in the alternating growth and decay 
of rolls that form the large scale circulation as discussed 
in [11. 



For Ra > 2.6 x 10^ all sdfs measured in the bulk have a 
wide range of active scales, and this range increases with 
Ra. These sdfs are continuous, and are not dominated 
by a small number of frequencies as one might expect in 
spatio-temporal chaos. Clearly these sdfs indicate that 
the bulk in the RB cell is turbulent for Ra > 2.6 x 10^. 
Although an inertial subrange is absent for the measure- 
ments at lower Ra, the sdf at Ra — 2.4 x 10^ (Fig. 6(j) ) 



shows an inertial subrange spanning roughly one decade. 
Note that the presence of an inertial subrange is not a 



7 




(i)Ra = 4.2 X 10** (LDA) (j)Ra = 2.4 x 10^ (LDA) 



Figure [S] (cont.) 



necessary condition for turbulence. An inertial subrange 
can only form when there is a range of scales which is 
neither influenced by macro- or micro-scales, i.e. at high 
Re when the production and dissipation scales are suf- 
ficiently separated. Therefore, an inertial subr ange i s 
absent despite the turbulence at lower Re [see also [g. [25j . 

The most striking observation is that the sdfs in the 
bulk and in the boundary layer are practically identical 
for Ra > 2.6 x 10^. As the sdfs for the bulk indicate that 
the flow is turbulent, one cannot escape from the con- 
clusion that the flow in the boundary layers must also 
be characterised as turbulent. It is difficult to reconcile 
this observation with the conceptual image of a laminar, 
Blasius-type boundary layer along the walls of the RB 
cell. For a Blasius boundary layer, the laminar-turbulent 
transition occurs at a relatively high Reynolds number, 
i.e.. Re* « 10^ — 10® depending on the intensity and 
nature of disturbances in the flow [26]. Here, Re* is a 
Reynolds number based on plate length and outer veloc- 
ity. However, even at the highest Ra considered in this 
study, the value of Re* based on the cell width and the 
wind velocity is only 3 x 10^, which is far below the crit- 
ical value for laminar-turbulent transition in a flat plate 
boundary layer. Thus, the turbulence in the boundary 
layer cannot be due to a shear-triggered transition (oc- 



curing inside the boundary layer) , but more likely due to 
the impingement and detachment of plumes (travelling 
into and out of the boundary layer), which is corrobo- 
rated by the great similarity between the bulk and the 
boundary-layer spectra. 

Naturally, one can always envision that the boundary 
layer only quantitatively behaves as a laminar Blasius 
boundary layer: the plumes introduce a time-dependence 
on time-scales that prevent pure laminarity. The under- 
lying assumption is that on the average, the thermals and 
plumes are passive when it comes to scaling of various 
important integral parameters of the kinematic bound- 
ary layer such as the friction factor Cf and the boundary 

1/9 

layer thickness A„. Indeed, Cf scales as C/ oc Re ' , 
as one would expect for a laminar boundary layer past a 
flat plate at moderate Ra [i,[ll[ni|23[. 

The scaling of the boundary layer thickness A„ is 
slightly more complicated. Assuming a laminar bound- 
ary layer, scales as A„ oc Re"^''^ [1^, which is related 
to Ra as A^ (x Ra~^^^ upon assuming Re oc Ra^^^. Ex- 
periments show that a scaling exponent —0.25 holds well 
for the boundary layer thickness of the sidewalls [1^ , but 
not for the scaling of the bottom and top boundary layers 
0, [13, which scales with exponents from —0.11 to 
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(e)Ra = 1.6 X lO'' (LDA) (f)Ra = 2.6 x 10^ (DNS) 

FIG. 6: Spectral density functions for the horizontal velocity component u for points in the boundary layer (bl) and in the bulk 
of the RB cell. Ra ranging from 7.0 x 10* to 2.4 x 10^ The source of the data, DNS or LDA, is indicated. 



—0.16. The simulations in this paper confirm the weak 
Ra number dependence of (Fig. [2]). This disparity 
indicates that plumes impinging on and emerging from 
the boundary layer are not entirely passive, but actively 
influence the boundary layer thickness. Given the im- 
portant role of A„ in the theory of it seems to us 
that a detailed study of how the interplay of plumes, the 
wind and the boundary layer influences A„ is desirable. 
A first step in this direction is taken in a forthcoming 
study [H, [l^] , of which the results indicate that the scal- 



ing of All cannot be derived from laminar boundary layer 
theory and scales as A„ oc Ag^Re"^. 

Although the sdfs point towards the presence of tur- 
bulence in the boundary layer, the characteristics of the 
boundary layer are of a rather different nature than those 
of classical (i.e. shear-driven) turbulent boundary layers. 
Indeed, viscous forces are still significant in the bound- 
ary layer, as can be deduced from the domain height H 

in viscous units 5y = vjur^ where u-y = ^ ^fjlw is the 

friction velocity [2^. Assuming that ~ U/Xu-, the 
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(i)Ra = 4.2 X 10* (LDA) (j)Ra = 2.4 x lO'* (LDA) 



Figure |6] (cont.) 



ratio H/5v which is the shear-Reynolds number Rer , can 
be approximated by 



(3) 



The simulations show that Re^. — 70 at Ra = 2.1 x 10^, 
and the viscous thickness of the boundary layer Xu/Sy — 
Xu/HKCr is only a fraction of this already very low value. 
Therefore, viscous effects are significant in the boundary 
layer. 

It is interesting to note that the scaling of the bound- 
ary layer thickness for flat plate flow changes from S/x (x 
Re~"^/^ (laminar) to S/x cx Re"^'''^ (turbulent, see [soj). 
where Re^ is the Reynolds number based on the distance 
X to the leading edge of the plate, and the thickness 6 
is the location where the velocity is 99% of the outer ve- 
locity. The decreased exponent resembles larger entrain- 
ment and thus a quicker growth of the boundary layer as 
a function of x. Upon assuming that the typical develop- 

1 /r) 

ment length is H and that Re oc Ra ' , the anticipated 
scaling for a turbulent boundary layer is S/H (x Ra~^^^^. 
If a « —0.11 is the scaling exponent A„ oc Ra" from the 
simulations, then — 1/4 < a < —1/14, so the exponent 
may resemble an intermediate state. 



However, it is not clear whether the entrainment mech- 
anism of a developing boundary layer past a flat plate is 
transferable to RB convection; the active role of plumes 
which transfer fluid (and heat) into and out of the bound- 
ary layer may significantly affect the entrainment charac- 
teristics. It is our suspicion that the relatively weak forc- 
ing in the horizontal direction due to the wind (which is 
responsible for the friction factor C/) is only weakly cou- 
pled to the relatively strong forcing in the wall- normal di- 
rection due to buoyancy. This may be the reason why the 
boundary layer has an sdf which is typical of developed 
turbulence but has many laminar features at the same 
time. Due to the weak coupling, the conceptual image 
of a laminar boundary layer may be appropriate for Cf 
but not for A„ in the boundary layer in RB convection. 
Therefore, the straightforward application of well-known 
characteristics of the Blasius-type boundary layer (such 
as the criteria for laminar-turbulent transition and the 
scaling of A„ with Re) to the boundary layers in RB con- 
vection is questionable without additional justification. 
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IV. CONCLUSION 

A combined numerical and experimental investigation 
was carried out to study the boundary layers in RB con- 
vection in a 4:1 aspect-ratio cell filled with water. The 
study covers a wide range of Ra numbers varying be- 
tween 7.0 X 10^ and 2.4 x 10^. The results for the lower 
Ra numbers were obtained from direct numerical simu- 
lations, while those for the higher Ra numbers followed 
from laser Doppler anemometry. In the midrange both 
techniques have an overlap of about one decade. 

The probability density functions of the horizontal ve- 
locity components have a Gaussian-like shape for the 
complete Ra number range, except at the lowest Ra 
where the flow is in a state of spatio-temporal chaos. 



Probability density functions that are scaled with the ve- 
locity rms value overlap for the entire Ra number range, 
with the exception of the largest Ra numbers where the 
large-scale circulation, or "wind" , is predominantly visi- 
ble in the histograms. 

The spectral density functions of the horizontal veloc- 
ity component that were determined in the bulk flow and 
in the boundary layer are surprisingly similar. Spectra in 
both regions of the flow are continuous and have a wide 
range of active scales for Ra > 2.6 x 10^. This indicates 
that both the bulk and the boundary layers are turbulent 
in the Ra number range considered, while (i) molecular 
effects are still noticeable in the boundary layer and (ii) 
the boundary layer does not behave like a classical shear- 
driven turbulent boundary layer. 
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